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Abstract 


We perform Monte Carlo simulations of a gauge invariant spin system 
which describes random surfaces with gonihedric action in four dimensions. 
The Hamiltonian is a mixture of one-plaquette and additional two- and three- 
plaquette interaction terms with specially adjusted coupling constants |I|, 

For the system with the large self-intersection coupling constant k we observe 
the second-order phase transition at temperature Pc — 1-75. The string ten¬ 
sion is generated by quantum fluctuations as it was expected theoretically . 
This result suggests the existence of a noncritical string in four dimensions. 
For smaller values of k the system undergoes the first order phase transition 
and for k close to zero exhibits a smooth crossover. 


1 Introduction 


The gonihedric string has been dehned as a model of random surfaces with an action 
which is proportional to the linear size of the surface |^, ||, 


A{M) = ^ Xij ■ Q{aij), 0(a) = Ivr - (1) 

<ij> 

where Ajj is the length of the edge < ij > of the triangulated surface M 2 and a^- is 
the dihedral angle between two neighbouring triangles of M sharing a common edge 
< ij >. The angular factor 0 dehnes the rigidity of random surfaces and for <^ < 1 
the angular factor increases sufficiently fast near angles a = vr to suppress transverse 
fluctuations [Q . The gonihedric action has been dehned for self-intersecting surfaces 
as well 0. The action accounts self-intersections of different orders by ascribing 
weights to self-intersections. These weights are proportional to the length of the 
intersection multiplied by the angular factor which is equal to the sum of all angular 
factors corresponding to dihedral angles in the intersection P|.The coupling constant 
in front of this term is called self-intersection coupling constant |^. In principle 
this coupling constant is a free parameter of the theory, but if one applies the 
continuity principle then one can £x the value of fee = 1/2 [y, Q (see formulas (4) 
and (28) in P|). 

The model has a number of properties which bring it very close to the Feynman 
path integral for a point-like relativistic particle P, This can be seen from (|I|) 
in the limit when the surface degenerates into a single world line, and in that case 
the action is proportional to the length of the path. This property of the gonihedric 
action guarantees that the spike instability, which is common to other triangulated 
random surface theories, does not appear here, because the action is proportional 
to the total length of the spikes and thus suppresses the corresponding huctuations 

The other important property of the theory is that at the classical level the 
string tension is equal to zero and quarks viewed as open ends of the surface are 
propagating freely without interaction 0 . This is because the gonihedric action (|l|) 
is equal to the perimeter of the hat Wilson loop |p]. As it was demonstrated in p, 
quantum huctuations generate a nonzero string tension 


(^quantum 2 

CL L) 


( 2 ) 


where d is the dimension of the spacetime, (3 is the coupling constant, a is a scaling 
parameter and 


d-2 

d 

in (p. In the scaling limit (3 ^ (dc = dje the string tension has a finite limit while 
the scaling parameter tends to zero as 


a = {l3- I3cf/^, 



thus the critical exponent v is equal to one half 


z/ = XjdH = 1 / 2 , 


( 3 ) 


where du is the Hausdorff dimension. Thus, although at the tree level the the¬ 
ory describes free quarks with string tension equal to zero, quantum fluctuations 
generate nonzero string tension and, as a result, the quark conhnement [^. The 
gonihedric string may consistently describe asymptotic freedom and conhnement, as 
it is expected to be the case in QCD[]. 

In addition to the formulation of the theory in the continium space the 
system allows an equivalent representation on Euclidean lattice Z'^ where a surface 
is associated with a collection of plaquettes. Lattice spin systems whose interface 
energy coinsides with the action (P have been constructed in arbitrary dimension 
d [|^ for the self-intersection coupling constant k = 1 and for an arbitrary fc in |^. 
This gives an opportunity for numerical simulations of the corresponding statistical 
systems in a way which is similar to the Monte Carlo simulations of QCD []^ . 

In three dimensions the corresponding Hamiltonian is equal to |l|, 


gonihedric 


= —2k ^ + — X/ ^rO'r+a+P o ^ ('^) 


2 Z_/ " ' '^r+a+P 2 

r,a,/3 r,a,l3 


and is a natural extension of the 3D Ising model 

Lj3d _ _ ^ ^ 

^Ising~ / . '^r'dr+a 


to the gonihedric case. The extensions of the 3D Ising model have been considered 
in the literature and their phase structure has been investigated. The essential 
point is that the geometrical nature of the gonihedric system specihes the coupling 
constants and the symmetry of the system^. This is in analogy |l|, ^ with super- 
symmetric systems where the symmetry also specihes the coupling constants. 

The self-intersection coupling constant k dehnes the degree of degeneracy of the 
vacuum state [^, |^; if /c 7 ^ 0, the degeneracy of the vacuum state is equal to 3 ■ 2^ 
for the lattice of the size and is equal to 2^^ when k = 0. The last case is a sort of 
supersymmetric point in the space of gonihedric Hamiltonians. What is amazing is 
that not only the ground state, but all of the energy levels have the same exponential 
degeneracy. This enhanced symmetry at the point k = 0 allows to construct the 
dual Hamiltonian in three dimensions jlSl , [16[| . For k ^ 0 the degeneracy of energy 


levels is smaller compared with the vacuum state. This exponential degeneracy is 
reminiscent of spin glass systems. 


^ Different modifications of the action have been proposed in the literature by adding the 
area term or the gaussian term. In the first case this will cause a nonzero string tension already 
at the classical level and will bring back all problems of the critical string when applied to stroM 
interactions. In the second case the model becomes simply equivalent to the gaussian model 0 
(see formulas (46)-(48)). 

^Spoiling of this fine tuning of coupling constants is equivalent to an addition of an area term 
and thus leads to nonzero string tension at the tree level (see previous footnote) and ’’drives” the 
system to an unwanted class of universality. 
























Despite the fact that the system has a rich symmetry, in three dimensions the 
two-point correlation function remains as a fundamental observable 
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C{f) =< (To- (T^? > 

and its behaviour at large distances is the best indicator of the phase transitions in 
the 3D gonihedric spin system. As usual, the most direct indication of the second- 
order phase transition is the growth of the correlation length, which in the scaling 
limit should tend to inhnity. The singularities of the energy, of the specihc heat 
and of the generalized magnetization are also indicative of the order of the phase 
transition. The generalized magnetization is dehned as the vacuum expectation 
value of the projection of spin states to the one which corresponds to the vacuum 
spin conhguration 


17 


=< (7^(vac) 


■ > 


(5) 


where a^{vac) denotes the vacuum spin conhgurations and /i = 1,2, for the 

case k = 0 and p = 1, 2,..., 3 ■ 2^ in all other cases. The exponential degeneracy of 
the ground state is common to gonihedric spin systems and the analogy with spin 
glasses allows to construct the generalized magnetization (^) . In the real exper¬ 

iments one can choose some subset of generalized magnetizations and measure 
the projection of virtual states to that subset of order parameters . 


The hrst Monte Carlo simulations 0 demonstrated that the gonihedric system 
with intersection coupling constant equal to one, k = 1, undergoes the second- 
order phase transition at f3c ~ 0.44. This happens near the critical temperature 
of the two-dimensional Ising model (3c = ^ln{l + ~ 0.44 |^. This result 

follows from the transfer matrix approach [T^ which describes the propagation of 
the two-dimensional system with the length and curvature amplitude. The system 
has a continuum limit at this temperature and describes two-dimensional free Dirac 
fermion [^. At the same time it was not expected that the critical indices will 


coincide with the ones of the 2D Ising model simply because the physical picture of 
the fermionic string propagation which follows from the transfer matrix approach 
is different from the propagation of the free point like fermion in two dimensions. 
In addition the analytical result (|) predicts the value z/ = 0.5, which does not 
coinside with the critical exponent of the 2D Ising model u = 1. The Monte Carlo 
simulations conhrm that the value of u is small, u ~ 0.34, and is almost twice 
smaller than that for the 3D Ising model where u = 0.63 In |]^ the value 
u ^ 0.44 is obtained. Thus the correlation function grows slower near the critical 
point. This means that the 3D gonihedric system and the 3D Ising model are in 
different classes of universality. 

Our aim is to study the phase structure of the spin systems which simulate 
random surfaces with gonihedric action in four dimensions. As it is well known, 
random surfaces with an area action in four dimensions can be simulated by the one- 
plaquette gauge invariant action (0) i- It should be reminded also that this gauge 
invariant spin system in four dimensions is self-dual and the critical temperature is 
equal to (3c = \ln{(i -|- \/2) |Q]. The Monte Carlo simulation of the system strongly 




















































indicates that the phase transition in 4D gauge invariant spin system is of the 
hrst order 

In the next section we shall describe the corresponding gonihedric Hamiltonian 
from [|l], H and the appropriate observables. The results of the Monte Carlo simula¬ 
tions are presented in the fourth section. For the system with the self-intersection 
coupling constant /c > = 1/2 we observe the second-order phase transition at 

temperature /dc — 1.75. The string tension is generated by quantum fluctuations as 
it was expected theoretically 0. This result suggests the existance of a noncritical 
string theory in four dimensions. 

For smaller values of k the system undergoes a hrst-order phase transition at 
temprature /3c > 2. The passage from second-order phase transition to the hrst-order 
one is caused by the condensation of self-intersections. Finally, for k < = 1/6 

the system exhibits a smooth crossover. Thus two critical values of k separate three 
regions of different critical behaviour and demonstrate the essential dependence of 
the phase structure on the intersection coupling constant k. 

2 Gonihedric spin system in four dimensions 

The system of random surfaces with gonihedric action can be formulated not only 
in three dimensions, but also in any dimension Similar to the three-dimensional 
case , the surface can be associated with a collection of plaquettes on a lattice and 
the interaction between spins can be organized in a way that the surfaces of interface 
will have the gonihedric action [^. 

In four dimensions the gonihedric Hamiltonian essentially differs from the three- 
dimensional one (ID because now it represents a gauge invariant spin system. The 
gauge invariant Hamiltonian in four dimensions has the form [Q 

HgLhedric = “4 II + H (aaaaaY^aaaaa), (6) 

^ {plaquettes} ^ {right angle plaquettes} 

where g‘^ is the gauge coupling constant and the independent spin variables aij should 
be attached to the centers of the edges < i,j > of the four-dimensional lattice. 

The Hamiltonian which simulates random surfaces with area action in four di¬ 
mensions is well known [Q and represents a gauge invariant spin system with one- 
plaquette interaction term 
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Haia = “4 H {aaaa). (7) 

^ {plaquettes} 

Thus the gonihedric Hamiltonian (|D in addition to the ’’ferromagnetic” interaction 
(1^ on each elementary plaquette has the ’’antiferromagnetic” interaction between 
two plaquettes which form a right angle [|^ 


J, 


plaquettes 


= -16 ■ J, 


rt plaquettes 


The partition function corresponding to (^) is equal to 





■0 g^H/A 


zm = i: e- 

{^r} 

These formulas completely define the gonihedric system in four dimensions and allow 
to simulate random surfaces on a four-dimensional lattice. It is important to stress 
that both systems and ( 0 ) are of a geometrical nature because in the hrst case 
the amplitudes are proportional to the linear size of a surface and in the second case 
they are proportional to the area. The other difference between these two systems is 
that in the gonihedric case (j^) the self-intersections of the surface can be propertly 
counted P, There are various formulations of self-avoiding random surfaces on 
a lattice |p. The main difference between these models and the gonihedric case is 
that the energy ascribed to self-intersections essentially depends on the geometry of 
intersecting plaquettes. 

Indeed, the essential property of the spin realization of the gonihedric system (P) 
is that it properly counts the self-intersections On a four-dimensional lattice a 
two-dimensional closed surface can have self-intersections of different orders because 
at a given edge one can have intersections of four or six plaquettes [|I|. The energy 
ascribed to a self-intersection essentially depends on a configuration of plaquettes 
in the intersection and there are two topologically different configurations of 
plaquettes with four intersecting plaquettes and only one with six intersecting pla¬ 
quettes. The corresponding energies as they are defined by the Hamiltonian (P) 
= 1) are equal to [Q 

62 = 1, 64 = 4, 64 = 5, 66 = 12 (8) 

and are equal to the number of plaquettes which intersect at a right angle. The 
total energy of the surface which has n 2 edges with two intersecting plaquettes at a 
right angle - ’’simple edges’’^, ^4 and h 4 edges with intersection of four plaquettes 
and hq with six plaquettes is equal to 

6 = n2 + 4h4 -|- 5fi4 -|- 12?7,6 (9) 

In the article |0] the intersection coupling constant k was introduced to control the 
intensity of self-intersections. In that case the Hamiltonian has the form 

51 hh 

^gonihedric 2 ^ ^ ^ ^ ^ {o'dCTCTa) 

^ {plaquettes} ^ {right angle plaquettes} 

- {aaaaaY^cTacraagY^agaaa). ( 10 ) 

^ {triples of right angle plaquettes} 

As it can be seen from the last expression in contrast to the case k = 1 (I) there 
appears a new three-plaquette interaction term^ 

The corresponding energies ascribed to self-intersections are equal to [Q 

■^edges with two intersecting plaquettes at a right angle 
The coefficient in front of this term is twice smaller than that in because for convinience 
we take all symmetric combinations of the plaquette triples inside the 3d cube. 





e2 = 1, €4 = 4fc, 64 = 6k — 1, eg = 12k 

and coinside with the weights (P) when k = 1. The total energy of the snrface in 
this general case is 

e = n 2 + Akn^ + {6k — l)n 4 + 12knQ ( 11 ) 

and rednces to when k = 1. 

Unfortnnately the method used for the construction of the transfer matrix in 
three dimensions 0 can not be used in four dimensions and we face an even a more 
complicated problem. 

3 Observables 

In pure gauge theories the magnetization is always equal to zero, < a >= 0, therefore 
this quantity can not be the order parameter. In other words the system does 
not have local order parameters to distinguish the phases and only gauge invariant 
products of spin variables over closed loops can be nonzero 

w{c) = < > ( 12 ) 

c 

and, as it was pointed out by Wegner 0, can separate the phases of the system. 
The perimeter and the area law of (|T^ can distinguish the ordered and disordered 
phases [|] (see formulas (3.34) - (3.39) in [Q). The same conclusion is correct for 
the gonihedric systems (^ and (0), therefore only different loop products are to be 
considered. 

To dehne the phase structure of the four-dimensional gonihedric gauge system 
one can measure different observables. The simplest one is the total energy per edge 
which can be expressed as a combination of different loop products of spins inside a 
three-dimensional cube [Q 


_ 1 _ LT 

Hidge = 3k - —— ^ {aaaa) + ^ {(raaa^Y\ao,aaa) (13) 

{six} {twelve} 

~ ^ {(^(^(^(^aT\craCraa/3Y\ai3aaa). 

{twelve} 

Here the hrst term denotes the product of spins around one plaquette in a 3d cube, 
the second one - a product of spins around two perpendicular plaquettes and the 
last one - a product of spins around three perpendicular plaquettes inside a 3d cube 
of the four-dimensional lattice. 

The advantage of this formulation of the theory is that one can measure the part 
of the surface energy which corresponds only to a simple edge which is proportional 
to n 2 

= 6 + ^(aaaa) - 2 ^ (aaaa)"(aaaa) (14) 

{three} 
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Htotal, 10'*4, fixed boundaries, kappa=1.00 


Figure 1: The average energy per edge (0) as a function of /3 . The system was 
heated and then cooled on steps of 0.001, points are plotted every fortieth step, 
k = l. 


- {aaaaaY\craCraai3y\apaaa). 

{twelve} 

The separation of the simple edge energy allows to compute the part of the surface 
energy which corresponds to self-intersections 

Itself—inter sections Tjtot Tjsimple 

^edge ^edge ^edge 

One can express this observables directly in terms of numbers of intersections 


Tjtot 

edge 


77-2 -|- 4/c774 -|- {6k — 1)?t.4 -|- 12knQ 


Tjsel f—inter sections 
edge 


rr simple 
^edge ^2 

4^774 -|- {6k — 1)77,4 + l2knQ 


The other important observable is the surface area which can be computed by using 
one-plaquette product 


Fedge 3 ^ ^ ) ((JCTCTCt), 


(15) 


where {aaaa) = P. The advantage to work with these observables is that they allow 
to control the inportance of the self-intersections in our attempt to hnd a nontrivial 
theory in four dimensions. 


4 Monte Carlo simulations 


As it is well known, random surfaces with an area action can be simulated by the 
gauge invariant action (|^) [Q. It should be reminded that this gauge invariant spin 
system in four dimensions is self-dual and the critical temperature is equal to [3c = 
4 / 77(1 -|- ^/2) 1^. The Monte Carlo simulation of the system strongly indicates that 
the phase transition in 4D Z 2 gauge invariant spin system is of the hrst order [0, n . 
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Htotal, 10**4, fixed boundaries, kappa=1.40 


Figure 2: Total energy on edge (p!3D , k = 1.4 
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Figure 3: Simple edges k = 1.4 



Figure 4: Self-intersection edges, A; = 1.4 












Figure 5: Plaquette, k = 1.4 


This has been indicated by measuring one-plaquette average (h5|) in the thermal 
cycle. A thermal cycle of the statistical system provides a general overview of its 
phase structure and is very helpful in defining the regions of the phase transitions 
[|]. Clear hysteresis in 4D Z 2 gauge invariant spin system is strongly indicative of 
the first-order phase transition [Q. The hysteresis is due to the metastability of the 
ordered phase at high temperature and of the disordered phase at low temperature. 

The similar Monte Carlo simulations of the gonihedric system (j^) and (Eg) on 
the lattices of the sizes 4"^, 6^, 8^, and 10^ demonstrate that the system has three 
regions of different critical behaviour which are characterized by the value of k. The 
thermal cycles with steps 5/9 = 0.001 have been performed in a large interval of the 
self-intersection coupling constant k from zero to four and we see that the phase 
structure essentially depends on k. 

For the system with the self-intersection coupling constant k > k^ = 1/2 we 
observe the second-order phase transition at temperature /9c — 1.75 (see Figures 
1, 2). figures 1, 2 show the energy per edge as a function of (3 in Monte Carlo 
simulations. The energy is a continuous function of temperature for large values of 
the self-intersection coupling constant k. The critical temperature does not actually 
depend on k. The small hysteresis loop still seen on the figures is a remnant of the 
loop which shrinks when we change 5(3 from 0.1 to 0.001 and pass to large lattices, 
this is indicative for the second-order phase transition. The behaviour of the simple 
edges, self-intersections and plaquette averages are shown respectively on figures 3,4 
and 5. From the figures 3 and 4 one can conclude that only at high temperature 
/9 < 0.5 the energy is concentrated on the self-intersections and that at the phase 
transition point the amount of energy on the self-intersections is negligible. The 
behaviour of the plaquette average on figure.5 shows that the area of the random 
surfaces scales at the critical point and signals that the string tension is generated 
by quantum fluctuations, as it is expected theoretically (see (||)). The Creutz ratio 
Q will be presented in a separate place. 

For smaller values of k the system undergoes a first order phase transition at 
temeprature (3c >2 (see figures 6, 7 and 8 ). The energy is discontinuous function 
of temperature, see figure 6. The critical temperature depends on k and is equal 
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Hbend, 10’'4, fixed boundaries, kappa=0.30 


Figure 7: Hysteresis loop for simple edges, k = 0.3 
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Figure 8: Hysteresis loop for self-intersection edges, k = 0.3 
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Figure 9: Hysteresis loop for plaquette, k = 0.3 
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Hbend, 10’'4, fixed boundaries, kappa=0.10 

Figure 10: Simple edges, k = 0.1 

to Pc — 2.5 for k = 0.3. At the phase transition point the amount of energy on 
self-intersections considerably increases ( compare hgures 2, 3, 4 and hgures 6, 7, 8) 
and is comparable with the energy on simple edges ( see hgures 7 and 8). Thus the 
passage from the second-order phase transition to the hrst-order one is caused by 
the condensation of self-intersections. The reason why = 1/2 is a critical valne 
for the self-intersection conpling constant is that at this valne the self-intersection 
energy can be eqnivalently considered as being a snm of simple edges. The self¬ 
intersection dissolves into a number of simple edges and there is no snppression any 
more of self-intersections compared with the simple edges. The plaqnette average 
shows strong discontinnity, see Fig.9, and the area does not scale in this case. 

Finally, for k < kc =1/6 the system exhibits a smooth crossover, see Fig.10, 
and snrfaces ’’evaporate” at all temperatures. 







5 Conclusion 


In conclusion we would like to stress that for the system with the large self-intersection 
coupling constant k we observe the second-order phase transition at temperature 
/3c — 1.75. The string tension is generated by quantum fluctuations as it was ex¬ 
pected theoretically [^. This result suggests the existence of a noncritical string 
theory in four dimensions. 
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